2 Stability of solutions to ODEs 



How can we address the question of stability in general? 

We proceed from the example of the pendulum equation. We reduce this 
second order ODE, 



+ ysin0 = O, 

V 



to two first order ODE's. 



Write x\ = 6, x 2 = 0. Then 

X\ = X2 

9 . 

X2 = — -smxi 

The equilibrium points, or fixed points, are where the trajectories in phase 
space stop, i.e. where 

*-(S)- ff 

For the pendulum, this requires 

X2 = 

X\ = ±77,7T, n = 0,1,2,... 

Since sinxi is periodic, the only distinct fixed points are 

(5)-C) - (!)-(:) 

Intuitively, the first is stable and the second is not. 
How may we be more precise? 



2.1 Linear systems 

Consider the problem in general. First, assume that we have the linear system 

u\ = a n ui + CI12U2 

U2 = CL2\U\ + CL22U2 
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or 

u = Au 

with 

u(t)= ( Ul{t) \ and A=( an au ] 

Assume A has an inverse and that its eigenvalues are distinct. 
Then the only fixed point (where u = 0) is u = 0. 

The solution, in general, is 

= aie Xlt c{ + a 2 e A2i cl 

where 

• Ai, A2 are eigenvalues of A. 

• c[, 62 are eigenvectors of A. 

• a\ and (X2 are constants (deriving from initial conditions). 
What are the possibilities for stability? 

1. Ai and A2 are both real. 

(a) If Ai < and A2 < 0, then u(t) — > as t — > 00. 
^> stable. 

u 2 
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(b) If Ai > and A2 > 0, then u{t) — > 00 as t — > 00. 
=4> unstable. 

u 2 




(c) If Xi < < A 2 , 

• If u(0) is a multiple of c[, then u(t) — > as t — > 00. 

• If u(0) is a multiple of c~2, then u(£) — > 00 as t — > 00. 
=4> unstable saddle. 




2. Ai, A2 are both complex. Then 

A = a ± zg. 

Assuming is real, 

= e CTt (/3i cos ^ + sin qt) 

(Pi, 02 are formed from a linear combination of of A's eigenvectors and the initial conditions). 

There are three possibilities: 

(a) Re{A} = a > =4> unstable. 



m 2 
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(b) a < => stable. 




u 



u 



l 



(c) a = =^> marginally stable. 



'2 




7 



We leave the case of repeated eigenvalues to Strogatz (pp. 135-6). 
2.2 Nonlinear systems 

We are interested in the qualitative behavior of systems like 

x\ = fi(xi,x 2 ) 

where f\ and f 2 are nonlinear functions of x\ and x 2 ■ 



Define U{ = X{ — x* to be a small departure from the fixed point. 
Perform a Taylor expansion around the fixed point. 

For a one dimensional function g(x) we would have 




g(x* + u) ~ g{x*) + g'(x*) • u 
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Here we obtain 



fi(x h X 2 ) = fj {xl,X* 2 \ + T^(X*,XZ)UI + T^(xl,xl)u 2 + 0{u 2 ) 
=0 since fixed pt 

The first term vanishes since it is evaluated at the fixed point. 



Also, since 
we have 
Substituting ili 

where 



Ui = Xi = f i (x 1 ,x 2 ) 
fi(xi,x 2 ) above, we obtain 

u ~ Au 



A = 



( dh dh\ 

dx\ dx2 



x=x" 



VJ2 <JJ 2 

\ dx\ dx 2 / 

A is called the Jacob ian matrix of / at x*. 
We now apply these results to the pendulum. We have 

x\ = fi(x h x 2 ) = x 2 

X 2 = f2(xi,X 2 ) = 



g . 

■— sinii 

V 



and 



A = 



1 

-g/l 



for 



Jit 

x 2 








There is a different A for the case 



x 



x 



7T 





(The sign of g/l changes.) 



The question of stability is then addressed just as in the linear case, via 
calculation of the eigenvalues and eigenvectors. 
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